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ABSTRACT 



The differential equations and boundary conditions describing the 
periodic-flow heat exchanger (frequently referred to as rotary regen- 
erator) are sufficiently complicated to preclude complete analytical 
solution. Because of the interest in this type of heat exchanger con- 
siderable effort has been expended to obtain useful solutions to the 
analytical problem of predicting its performance. Several approximate 
solutions have been proposed but each of these is limited due to the 
restricted ranges of parameters in which they apply. 

A numerical finite -difference method of calculating the spatial 
temperature distribution in the periodic -flow type heat exchanger is 
presented by considering the metal "stream" in cross -flow with each of 
the fluid streams as two separate but dependent exchangers. In the 
development no assumptions are made which would restrict the range of 
parameters in which the analysis would be applicable. The exchanger 
effectiveness has been evaluated over the following ranges of dimen- 
sionless parameters: 

1*0 ^ ^min/ Cmax 0.10 

1.0 ^ Cr/Cmin 5? c* 3 

1.0 ^ NTU 0 ^ 10 

1.0 ^ (hA)' 0.25 

The National Cash Register 102A general purpose digital computer 
was used to accommodate the large number of subdivisions necessary for 
accuracy and an extrapolation of data to zero element area (infinite 
subdivision) was used to obtain values of effectiveness good to four 
significant figures. 

The author expresses his appreciation to Assistant Professor C.P. 
Howard for his direction and encouragement in this work. 
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NOMENCLATURE 



SYMBOLS 

A = heat transfer area on side designated by subscript, sq ft 
C = heat capacity rate (Wc) of fluids or rotor matrix according to 
subscript, Btu/hr, deg F 

c = specific heat of fluids (constant pressure) or rotor matrfx 
material depending on subscript, Btu/lb, deg F 
h = unit conductance for thermal convection heat transfer, 

Btu/hr, sq ft, deg F 

k = unit thermal conductivity, Btu/hr, sq ft, deg F/ft 
K = a constant, defined a3 used and distinguished by numerical 
subscript 

N » number of subdivisions of fluid or matrix solid flow "stream" 
according to subscript 
Q * heat transfer rate, Btu/hr 

* 

T = temperature of fluid or matrix, inlet or outlet, depending on 
subscripts, deg F 

W = mass flow rate of fluid (lb /hr) or matrix (rev /hr) (mass) 
according to subscript 
DIMENSIONLESS PARAMETERS 

E = exchanger heat transfer effectiveness; ratio of actual to 

thermodynamically- limited maximum possible heat transfer rate 
^min/Pmax = capacity rate ratio of fluid flow streams 
CrA-min = capacity rate ratio of rotor matrix to minimum fluid 
capacity rate 

(hA)' = (hA) n /(hA) x , geometric symmetry and heat transfer balance 

of exchanger 
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NTU = (hA)/C, number of transfer units on side designated by 



SUBSCRIPTS 

avg = average 
r = rotor matrix 
i = inlet 

o = outlet, or over-all in case of NHJ 0 
min = minimum magnitude 
max = maximum magnitude 
n = value on side of 
x = value on side of C^av 
c = cold side 
h = hot side 

° o = subscript on E to indicate value when Cr/Cmin = oO 



subscript 







over- all number of transfer units 
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1. Introduction 



In the past decade considerable interest has been shown in the ap- 
plication of the periodic-flow heat exchanger (rotary regenerator) to the 
Gas-turbine cycle for exhaust-gas thermal- energy regeneration. Regen- 
eration by whatever method applied to the gas-turbine cycle obtains, among 
other improvements, higher thermal efficiency, lower specific fuel con- 
sumption, improved part load performance and lower optimum pressure ratio. 
Particular advantages to be expected from the application of the periodic- 
flow heat exchanger to the gas-turbine process are discussed by Harper 
and Rohsenow (l) 1 and Coppage and London (2), among others. Principal 
among these is the possibility of a more compact unit than is possible 
with the conventional stationary surface types. This is due to the fact 
that the periodic-flow heat exchanger matrix can be finely divided, with- 
out structural difficulty, resulting in an appreciably reduced heat ex- 
changer volume relative to a stationary type for a given effectiveness 
and pressure drop. The major disadvantage of the rotating periodic-flow 
type is the sealing difficulties encountered and the attendant leakage 
involved when an appreciable pressure ratio is used. This has been re- 
duced to acceptable limits in at least two current development programs 
(3,4) and further refinements should follow. 

Because of the interest in this type of heat exchanger there has 
been considerable effort expended to obtain useful solutions to the 
analytical problem of predicting its performance. The differential 
equations and boundary conditions describing the system are sufficiently 
complicated to preclude complete analytical solution. An analytical 



lumbers in parentheses refer to the Bibliography at the end of the thesis. 
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solution is available for the special case of infinite matrix rotative 

speed vhere the expression for effectiveness takes the same form as that 

of a counter-flow direct type (2) and is given by 

-rrUo C i — C /C»V1 A X ) 

r / - C r j 

°° j _ C«\^/c mU)i 0 ” ( n " n /Cm6<) u j 

and for the further condition of C min/ C max * unity, reduces to 

V 



E 



OO 



A/rq p 
I +■ Nrtu 




Consequently, several approximate methods of solution have been 
proposed. Perhaps the oldest is that of Hausen (6), described and elab- 
orated by Johnson (8,9), which begins with the known solution for the 
initial cooling of a matrix element, initially at uniform temperature, 
by the passage of a cold fluid. By dividing the regenerator tube into a 
large number of strips of equal length and applying this elemental solu- 
tion to each strip in turn, a family of simultaneous linear algebraic 
equations may be derived whose solution gives the desired temperature 
variation in the periodic case. 

A second method, due to Iliffe (10), uses the solution to the prob- 
lem of the initial cooling of a matrix element at arbitrary initial temp- 
erature distribution, derived by Nusselt. Using this elemental solution, 
the method of Iliffe leads to two simultaneous integral equations which 
require numerical methods for their solution. In a third method, Saun- 
ders and Smoleniec (7) replace the governing partial differential equations 
by their equivalent finite difference equations end effect an approximate 
solution by relaxational techniques. In a fourth approach, Tipler (16) 
replaces the thermal problem by its electrical analogue end observes the 
solution on an oscilloscope. In yet a fifth approach, the perturbation 
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method used by Schultz (8) has been extended by Jones and Fax (5) lead- 
ing to algebraic statements for the departure of the regenerator effec- 



types of asymmetry. Each of these solutions is limited in one respect 
or another, primarily due to the restricted ranges of parameters in which 
they apply. 

It is the purpose of this thesis to present the results of a finite- 
difference numerical analysis which was carried out with the aid of a 
digitial computer such that there were no limitations imposed on the four 
independent dimensionless parameters selected. 



tiveness from its asymptotic value, given by equation 
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2. Method 



The conventional idealizations and boundary conditions assumed in 
the usual derivation of the governing differential equations (2) are as 
follows : 

a. The thermal conductivity of the matrix is zero in the direction 
of fluid flow and matrix metal "flow". It is infinite in the 
other direction normal to the fluid flow. 

b. The specific heats of the two fluids and matrix material are 
constant with temperature. 

c. No leakage of the fluids occurs either due to direct leakage 
or carry over, and each fluid flow is unmixed. 

d. The convective conductances between the fluids and the matrix 
are constant with flow length. 

e. The fluids pass in counterflow directions. 

f . Entering- fluid temperatures are uniform over the flow inlet 
cross section and constant with time. 

g. Regular periodic conditions are established for all matrix 
elements, i.e., steady state condition. 

With these idealizations a typical element, Fig. 1, from each fluid 
side of the exchanger can be represented schematically as shown in Fig. 

2. Although Fig. 1 is a disk (axial-flow) type exchanger, typical el- 
ements of a drum (radial-flow) type can be selected, with the same ass- 
umptions, which can be represented schematically in the same manner as 
Fig. 2. 

It is apparent from Fig. 2 that each of these elements can be regard- 
ed as a crossflow exchanger with a gas stream and a metal "stream". This 
extension of the method of calculating crossflow exchangers to the rotary 
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FIG. 2 SCHEMATIC HEAT EXCHANGE ELEMENTS 



regenerator was suggested by G. M. Dusinberre in his discussion of ref- 
erence ( 2 ). It should be noted that the T rl and T r0 on the two elements 
are not the same, but to avoid further subscript notation the side in 
question will be specified when it is not obvious. However, the T ro of 
a particular element is the T r ^ of the next adjacent element in the 
direction of matrix flow, and this is of particular significance at the 
seals where the direction of fluid flow reverses. 

Consider for the moment only the heat exchange element on the side 
of CniAy . Fig 2(a). For simple unmixed cross flow with uniform tempera- 
tures across the inlets and average temperatures at the outlets, the 



heat transfer rate for this element is 

Q = ^max(^xi " T xo )( 1/N x ) [2] 

Q * C r (T r0 - T ri )(l/N r ) [ 3 ] 

Q = (hA) x ZiTaYgCl/NjJljf) [^1 



For a small enough element the arithmetic-mean temperature differ- 
ence may be assumed valid, so that 

/^ T avg = (l/2)(T x ^ + T xo ) ~ (l/2)(T r £ + T ro ) .... [*5j 
By using equations [2] , [5] , [k] , and [ 5J the outlet temperatures « 
of the two streams may be solved for in terms of the inlet temperatures 



and put in the form (Appendix 1) 

Txo ■ T xl - (T x t - T rl ) [6a] 

T r o = ^ri + Ks(T x i ~ Tj-j) . . . [a>J 

By the same method the outlet temperatures of the element on the 
side of Cjuin may be solved for in terms of its inlet temperatures giving 

Tno • T nI ♦ Ka(T rl - T nl ) [7a] 

T ro » T rI - K 4 (T rl - T nI ) M 



By using the definitions of the parameters and proper algebraic 
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manipulation the foregoing constants may be put directly in terms of 



the dimensionless parameters and the number of subdivisions of the matrix 
and fluid streams, giving 

“ 1-1 



K, = Z 
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Now consider the schematic representation of the composite heat 
exchanger made up of elements as shown in Fig. 3. For illustration the 
fluid streams and the matrix stream have each been divided into three 
equal substreams to form the elements. The double lines represent the 
area of the seals. It should be noted that the left edge is physically 
the same as the right edge and therefore the matrix inlet temperature of 
a particular matrix substream on the left must be identical to the matrix 
outlet temperature of that substream on the right. This will be referred 
to later as the reversal condition. 

In Fig. 3 where individual element temperatures are designated with 
the same label it is not implied that they are equal numerically, except 



8 




9 



FIG. 3 SCHEMATIC COMPOSITE HEAT EXCHANGER 



where indicated as constant across the fluid inlets. This is done to 
avoid additional subscript notation. Also, it is not to be implied from 
Fig. 3 that there is necessarily any symmetry with respect to the areas 
involved . 

Now for calculation purposes if the fluid inlet temperatures are 
given some convenient values, such as one and zero, and a matrix inlet 
temperature distribution is assumed on the left edge, the remaining 
temperatures of each element can be calculated in the order indicated 
by the element numbers by repetitive use of equations [6] and then lj . 
If the temperature distribution assumed on the left edge was correct 
then it would be duplicated on the right, i.e., the reversal condition 
was fulfilled, and the problem would be solved for the particular set 
of parameters used. If, however, this is not the case, then the result- 
ing temperature distribution on the right is now used on the left and the 
procedure repeated until the reversal condition is met. 

After several iterations of this type the cyclic character of the 
calculations for each element and column on a side becomes apparent and 
the beauty of automatic computing equipment is fully appreciated. It 
can be seen that if, in the process of computing the outlet temperatures 
of element one, the T x0 is put where the T xl was and a control number 
increased by one so that the T r ^ of element twowill be used next the 
same sequence of operations can now be used to compute element two, and 
so on down the first column of elements. Also, if the T r0 of each ele- 
ment is put where its T r ^ was, then upon reaching the end of the column 
and "zeroing” the above control number the same cyclic operation can be 
used for the next column. A second control number is used to indicate 
when one side is finished and then a similar procedure is followed on 
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the other side. In this scheme it is only necessary to provide stor- 
age for the T r i and T ro of each row, which will be used for comparison 
purposes to determine if the reversal condition is met, and the Txo of 
each column on one side and the T no of each column on the other side, 
which will be used to make a heat balance and compute the effective- 
ness when the reversal condition is met. 

That the reversal condition eventually will be met by this iterative 
scheme can be shown (Appendix 1 ) to depend on the physical conditions of 
the problem and the second law of thermodynamics. This will be called 
the condition for convergence and for this problem, using equal sub- 
divisions of each stream, it may be stated in the following form: 

NTU 0 [l + l/(hA)'J Cmin/Cmax 4 2N [ 9 ] 

Accepting for the moment the idealizations and assumptions made, 
then the accuracy of the solution obtained, by the very nature of a 
finite-difference process, depends only on the number of subdivisions 
used. It is apparent from condition [9J that a greater number of sub- 
divisions will also enhance the convergence. 

A considerable number of subdivisions can be handled by automatic 

«* 

computing equipment to obtain a certain accuracy and aid convergence, 
but the number of subdivisions used must be arrived at by a compromise 
between the accuracy desired and time available. The easiest way to 
illustrate the accuracy obtained for a given degree of subdivision is 
to solve the problem for several values of subdivision and extrapolate 
the results obtained to zero element area. As an example, consider the 
following results for the indicated parameter values: 

^min/^max = NIU 0 = 6, ^r/^min = = 
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No. subdivisions 



Effectiveness 



8 


0.7816 


16 


0.7800 


52 


0.7796 



These results are plotted versus element area and extrapolated to 
zero In Fig. 3a where it Is found that the indicated effectiveness is 0.7795* 
This correction is not constant for all ranges of parameters, being 
largest for low values of C r /C min and high values of NTU 0 , the fore- 
going example being typical. The procedure used in this work was to 
calculate Tables 1 through 9 with l6 subdivisions, then recalculate 
Tables 1 and 9 with 32 subdivisions. Table 9 was calculated to show 
the effect of varying the parameter (hA)', and it also contains rep- 
resentative values from the other tables so that with these data all 
of the tabulated data were extrapolated to zero element area. 

The computer used has the equivalent of nine significant decimal 
figure accuracy; it was specified that the reversal condition be ful- 
filled to only five figures. Of more significance than this was the 
final requirement that before a solution was considerd as being obtained 
the heat balance error, which is a direct measure of error in the effec- 
tiveness introduced by accuracy of computation, be less than 0.024 per 
cent. 

The remainder of the computer program was essentially a system of 
"feeding" it the next set of parameters after each solution. Appendix 
2 describes the computer program and its various modifications in some 
detail. Several interesting features incorporated in it are: 

a. Up to 64 subdivisions could be accommodated without external 
memory storage simply by modifying certain control commands. 
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b. Final matrix temperature distributions were punched on cards 
for future use. 

c. A logical choice of initial estimate for matrix temperature 
distribution for each new problem was made from data available 
from work up to that point. 

d. Optional read out of all or selected temperatures was provided. 
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3. Comparison of Results 

The results obtained are reported in tabular (Tables 1-9) and 
graphical form (Figs. 4 - 11). Several general observations can be made 
from the data: 

a. As C ^ p/C^y decreases, for given values of the other parameters, 
the effectiveness increases. This is predicted by the limiting solution 
for C r /C min = oo , equation [lj , and vas pointed out in the illustrative 
problem of reference (2). This is illustrated in Fig. 12. 

b. By using NTU 0 as defined in the nomenclature, which contains 
(hA)' explicitly, the influence of (hA)' at a given NTU 0 is very small 
for Cmin/CmHY ^ 0.9 . This conclusion was also reached by Coppage and 
London (2) in their evaluation of the work by Hausen (6), Saunders and 
Smoleniec (7)> J.E. Johnson (8), and Iliffe (10). For C min/ C max ^ 0.9> 
however, the influence of (hA)' becomes increasingly pronounced, and the 
variation for (hA)' = unity and 0.25 may be as much as seven per cent for 
Cmln/C niflY = 0.1. See Table 9» 

c. The difference in effectiveness for C r /C min ^ 5 and the value 
given by the limiting solution for Cp/Cmin = 00 is very small. This 
may be seen from inspection of Tables 1 through 8 or Figs. 4 through 11. 

R.W. Johnson (ll) used much the same method employed here but had 

limited access to a computer and was therefore able to calculate only a 

relatively small number of the points reported here. Also due to limited 

* 

computer time, at most ten and in some cases only five subdivisions were 
used with no extrapolation to zero element area correction. Allowing for 
this the corresponding results are practically identical. 

Coppage and London (2) recommended the results of J.E. Johnson ( 9 ) 
which were obtained for C m i n /C max = unity. These results, in terms of 
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the parameters used here, were later incorporated in reference (12) with 
extrapolation 'to Cmin/Cmax = 0.7 . The variation in corresponding values 
of Table 1 of this work and Table 6, page 11-32, of reference (12) which 
are for C m in/Cinax = unity is less than one per cent. Only when C m1n /C niflY 
= 0.7 and Cp/C^n = unity, where the extrapolation method used was ad- 
mitted to be most uncertain, is the variation slightly more than two per 
cent. 



An expression for effectiveness in closed form may be desirable for 
certain uses such as computer analysis of plant performance. A relative- 
ly simple empirical formulation is suggested in reference (12) for 
c min/ c max ^ 0.7 and Cp/Cmm ^ 2 and is of the form 



E oc — E 






A E 

E oo 



*5 (Cr/Cmln)" 




where Kg = l/9 , m = -2, and E^ is determined from equation [lj . Such 
a form cannot be applicable for all values of parameters with single val- 
ues of Kg and m, but it agrees reasonably well with the results reported 

here for the ranges specified. It has been found that with slight var- 

« 

iation of Kg and m very close agreement can be obtained in a small range 
of parameters which may be of particular interest. As an example, in 
the range 



1.0 


>> 


^min/^max ^ 


0.90 


1.25 


<1 


C r / C min ^ 


5.0 


5.0 




NTU q ^ 


9.0 



with Kg = l/9 and m = -1.87, the error is less than one per cent. Fig. 
* 

13 illustrates the method used to determine Kg and m in this range. 
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k. Limitations 



The most questionable of the idealizations made is that of zero 
matrix conductivity, particularly in the direction of fluid flow. Al- 
though the effect of matrix conductivity could be included in the finite- 
difference equations the additional degrees of freedom would make adequate 
coverage of all the parameters extremely voluminous and time consuming. 

If for a particular exchanger design the matrix conductivity was consider- 
ed important then the problem could be solved individually on a computer 
with adequate capacity. A simpler alternative is to neglect the conduct- 
ion and correct the effectiveness for this effect. From theoretical con- 
siderations Professor A.L. London has shown (13) that a correction factor 
for conductivity in the direction of fluid flow for approximately equal 
fluid capacity rates, of the form 

- J£_ [uj 

E LC 

should result in a somewhat pessimistic prediction of the reduction in 
effectiveness. (A in this case is the matrix cross-sectional area avail- 
able for heat conduction in the direction of fluid flow on both sides, 
and L is the flow length). Schultz (l4) has solved for the effectiveness 
as a function of the conduction parameter in equation [llj for three cases. 
In the nomenclature used here, these three cases sure for the following 
conditions: 

a. Cmin = very large, k finite 

b. = C,^ , Cy/C^ = 00 , k finite, regenerator sub- 

divided into sections in direction of fluid flow 

c. c min = c max , ^r/ c min finite, k = °o , regenerator sub- 
divided into sections in direction of fluid flow. 
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Prom the curves in reference (lk) which allow direct comparison with 
the results presented here it is found that for 



the reduction in effectiveness is at most about one per cent. Addition- 
al work, not available at the time of writing, concerning the influence 
of conduction is reported in reference (15). 

The effect of leakage can be calculated separately resulting in a 
correction to the effectiveness without leakage. This has been shown to 
be small for as much ten per cent leakage (l). 



kA / 0.02 

LC 
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5. Summary and Conclusions 



The numerical procedure used in this work assumes no symmetry , 
either with respect to heat capacity rates of the fluids or area ratios 
of the two sides, and therefore is not limited in this respect. The 
ranges of parameters which have been covered effectively are 

1.0 ^ ^min/^max ^ 0*10 

1.0 ^ C r /C mln ^ 

1.0 ^ NTU 0 ^ 10 

1.0 > (hA) ’ ^ 0.25 

The results are tabulated at frequent enough intervals and over suffi- 
cient ranges of the parameters to include the gas-turbine regenerator 
application as well as other applications, such as intercooling, where 
the values of the parameters might be quite different. 

The error introduced by the finite-difference type of solution has 
been minimized by use of sufficient subdivision and further extrapolation 
to zero element area so that four significant figures in the results are 
justified. 

The effects of deviations from the assumed idealizations have not 
been included in the results. 

Even though the results obtained cannot be made to fit a completely 
general closed form, a relatively simple expression for effectiveness 
can be obtained which will give good agreement over a reasonable range 
of parameters. 
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TABLE 1 



PERIODIC -FLOW HEAT EXCHANGER EFFECTIVENESS 



^mln/^max = 1*0 
(hA)' = 1.0 





1 




- C r /Cmin • 








1 


NTO 0 


1 


1.25 


1.5 


2 


5 


5 


00 


0.5 


0.5221 




O.5282 


0.5504 




0.5529 


0.5555 


1.0 


0.4665 


0.4779 


0.4845 


0.4912 


0.4960 


0.4985 


0.5000 


1.5 


0.5478 




0.5757 


0.5861 




0.5978 


0.6000 


2.0 


0.6007 


0.6251 


0.6560 


0.6491 


0.6587 


0 . 6658 


0.6667 


2.5 


0.6585 




0.6792 


0.6941 




0.7110 


0.7142 


5.0 


0.6672 


0.6959 


0.7119 


0.7280 


0.7400 


0.7465 


0.7500 


5-5 


0.6900 




0.7577 


0.7546 




0.7759 


0.7778 


4.0 


0.7086 


0.7410 


0.7586 


0.7760 


0.7889 


0.7959 


0.8000 


4.5 


0.7242 




0.7760 


0.7957 




0.8159 


O.8181 


5.0 


0.7575 


0.7725 


0.7907 


0.8086 


0.8217 


0.8290 


0.8555 


5.5 


0.7491 




0.8054 


0.8215 




0.8417 


« 0.8461 


6.0 


0.7592 


0.7956 


0.8144 


0.8525 


O.8455 


0.8526 


0.8571 


6.5 


0.7882 




0.8241 


0.8419 




0.8621 


0.8667 


7.0 


0.7765 


0.8158 


0.8527 


0.8504 


0.8652 


0.8704 


0.8750 


7.5 


0.7855 




0.8404 


0.8579 




0.8777 


0.8825 


8.0 


0.7901 


0.8285 


0.8475 


0.8647 


0.8772 


0.8842 


O.8889 


8.5 


0.7961 




0.8556 


0.8708 




0.8900 


0.8947 


9-0 


0.8017 


0.8407 


0.8595 


0.8765 


0.8884 


0.8955 


0.9000 


9.5 


0.8068 




0.8646 


0.8814 




0.9002 


0.9047 


10.0 


O.8115 


0.8510- 


0.8694 


0.8860 


0.8978 


0.9045 


0.9091 
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TABLE 2 



PERIODIC-FLOW HEAT EXCHANGER EFFECTIVENESS 



^minA'max = ®*95 



(hA)’ 



NTU 0 


1 

1 


1.25 


- ^r/^min 

1.5 


1 


0.4714 


0.4833 


0.4901 


2 


0.6082 


0.6319 


0.6454 


3 


0.6760 


0.7064 


0.7234 


4 


0.7183 


0.7528 


0.7715 


5 


0.7477 


0.7849 


0.8047 


6 


0.7698 


0.8090 


0.8292 


7 


0.7872 


0.8277 


0.8481 


8 


0.8012 


0.8429 


0.8633 


9 


0.8129 


0.8554 


0.8757 


10 


O.8229 


0 . 8661 


0.8861 



= 1.0 



2 


3 


5 


00 


0.4970 


0.5021 


0.5047 


O.5063 


0.6592 


0.6693 


0.6747 


0.6778 


0.7406 


0.7533 


0.7600 


0.7640 


0.7903 


0.8039 


0.8113 


0.8158 


0.8240 


0.8379 


0.8455 


0.8503 


0.8484 


0.8624 


0.8700 


0.8750 


0.8672 


0.8808 


0.8884 


0.893^ 


0.8820 


0.8953 


0.9027 


0.9077 


0.8940 


0.9069 


0.9142 


0.9191 


0.9040 


0.9165 


0.9235 


0.9284 



21 



TABLE 3 



PERIODIC-FLOW HEAT EXCHANGER EFFECTIVENESS 



^min/^max ~ 0.90 



(hA)' 





1 




- c r / c m±! 


NTU 0 


1 


1.25 


1.5 


1 


0.4763 


0.4887 


0.4958 


2 


0.6156 


0.6405 


0.6547 


3 


0.6846 


0.7167 


0.73^7 


4 


0.7275 


0.7642 


0.7841 


5 


0.7573 


0.7970 


O.8181 


6 


0.7795 


0.8216 


0.8432 


7 


0.7970 


0.8407 


0.8625 


8 • 


O.8112 


0.8561 


0.8780 


9 


0.8230 


0.8689 


0.8906 


10 


0.8331 


0.8797 


0.9012 



= 1.0 



2 


3 


5 


00 


O.503O 


0.5082 


0.5110 


0.5126 


0.6693 


0.6800 


0.6856 


0.6889 


0.7530 


0.7664 . 


0.7735 


0.7777 


0.8040 


0.8185 


0.8263 


0.8310 


0.8385 


0.8534 


0.8615 


0.8664 


6.8636 


0.8784 


0.8865 


0.8916 


0.8827 


0.8972 


O.9051 


O.9102 


0.8978 


0.9117 


0.9195 


0.9246 


0.9099 

0.9199 


0.9234 

0.9329 


0.9309 

0.9401 


0.9359 

0.9450 
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TABLE 4 PERIODIC-FLOW HEAT EXCHANGER EFFECTIVENESS 



C tmn/ G max = °- 8 ° 
(hA) ' = 1.0 



• 


r~ 




c r/ c min 






1 


NTO 0 


i 


1.25 


1.5 


2 


3 


00 


1 


0.4861 


0.4994 


0.5071 


0.5149 


O.5207 


0.5254 


2 


0.6299 


0.6573 


0.6730 


0.6891 


0.7011 


0.7109 


3 


0.7007 


0.7364 


0.7565 


0.7769 


0.7919 


0.8043 


4 


0.7443 


0.7853 


0.8077 


0.8299 


0.8461 


0.8597 


5 


0.7743 


0.8190 


0.8427 


0.8655 


0,8818 


0.8957 


6 


0.7967 


0.8440 


0.8683 


0.8910 


O.9069 


0.9206 


7 


0.8140 


0.8634 


0.8879 


0.9101 


0.9253 


0.9386 


8 

\ 


0.8279 


0.8788 


0.9033 


0.9249 


0.9392 


0.9518 


9 


0.8394 


0.8914 


0.9157 


0.9365 


0.9500 


0.9619 


10 


0.8492 


0.9021 


0.9260 


0,9460 


0.9586 


0.9696 
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TABLE 5 PERIODIC -FLOW HEAT EXCHANGER EFFECTIVENESS 



Cjnin/^max = 0.70 

(bA) ' = 1.0 





r 






NTU 0 


i 


1.25 


1.5 


1 


0.4959 


0.5104 


0.5187 


2 


0.6457 


0.6737 


0.6910 


5 


0.7155 


0.75^ 


0.7770 


4 


0.7590 


0.8044 


0.8291 


5 


0.7887 


0.8582 


0.8643 


6 


0.8104 


0.8629 


0.8895 


7 


0.8272 


0.8818 


0.9085 


8 


0.8405 


0.8967 


0.9232 


9 


0.8514 


0.9087 


0.9548 


10 


0.8606 


0.9188 


0.9443 



2 


5 




1 

00 


O.5271 


0.5555 




0.5584 


0.7087 


0.7219 




0.7526 


0.7995 


O.816O 




0.8295 


0.8535 


0.8710 




0.8855 


0.8891 


0.9063 




0.9207 


0.9159 


0.9503 




0.9459 


0.9520 


0.9475 




0.9598 


0.9456 


0.9596 




0.9709 


0.9560 


0.9687 




0.9788 


0.9642 


0.9756 




0.9845 
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TABLE 6 PERIODIC -FLOW HEAT EXCHANGER EFFECTIVENESS 



^mln/^max = 0 . 50 

(hA) ’ = 1.0 





1 




AWn 


NTO 0 


1 


1.25 


1.5 


1 


0.5156 


O.5322 


0.5417 


2 


O.669O 


0.7045 


0.7249 


3 


0.7402 


0.7869 


0.8133 


4 


0.7817 


0.8353 


0.8645 


5 


0.8089 


0.8670 


0.8975 


6 


0.8284 


0.8894 


0.9202 


7 


0.8431 


0.9062 


0.9366 


8 


0.8546 


0.9191 


0.9489 


9 


0.8639 


0.9293 


0.9582 


10 


0.8716 


0.9376 


O.9655 



2 


3 


! 

00 


0.5515 


0.5588 


0.5647 


0.7461 


0.7618 


0.7746 


0.8399 


O.859O 


0.8744 


0.8928 


0.9122 


0.9274 


0.9256 


0.9436 


0.9572 


0.9469 


0.9629 


0.9745 


O.9615 


0.9752 


0.9847 


0.9717 


0.9833 


0.9908 


0.9789 


0.9886 


0.9944 


0.9842 


0.9921 


0.9966 



* 
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TABLE 7 PERIODIC-FLOW HEAT EXCHANGER EFFECTIVENESS 



C mln/ C max = 0.30 

(hA) 1 = 1.0 





1 






NTU 0 


1 


1.25 


1.5 


1 


0.5346 


0.5558 


0.5648 


2 


0.6902 


0.7315 


0.7556 


5 


0.7577 


0.8114 


0.8422 


4 


0.7951 


0.8558 


0.8895 


5 


0.8191 


0.8838 


O.9183 


6 


0.8360 


0.9031 


0.937^ 


7 


0.8487 


0.9171 


0.9506 


8 


0.8588 


0.9278 


0.9602 


9 


0.8671 


0.9363 


0.9673 


10 


0.8740 


0.9433 


0.9729 



2 


3 


00 


0.5761 


0.5846 


0.5915 


0.7804 


0.7988 


0.8136 


0.8729 


0.8945 


0.9110 


0.9214 


0.9420 


0.9566 


0.9490 


0.9670 


0.9787 


O.9658 


0.9808 


0.9895 


0.9764 


O.9885 


0.9948 


0.9835 


0.9930 


0.9974 


O.9881 


0.9955 


0.9987 


0.9913 


0.9971 


0.9994 
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TABLE 8 PERIODIC-FLOW HEAT EXCHANGER EFFECTIVENESS 



C mln/ C max 3 °.10 

(hA)' = 1.0 





1 




C r /Cmin 


NTUq 


1 


1.25 


1.5 


1 


0.5554 


0.5752 


0.5877 


2 


0.7070 


0.75^2 


0.7820 


5 


0.7680 


0.8279 


0.8630 


4 


0.8009 


O.8667 


0.9042 


5 


0.8223 


0.8908 


0.9283 


6 


0.8378 


0.9075 


0.9440 


7 


0.8498 


0.9200 


0.9550 


8 


0.8595 


0.9298 


0.9632 


9 


0.8675 


0.9376 


0.9694 


10 


0.8743 


0.9441 


0.9743 



2 


3 


00 


0.6008 


0.6l04 


O.6186 


0.8109 


O.8319 


0.8487 


O.8981 


0.9219 


0.9391 


0.9397 


0.9615 


0.9753 


0.9618 


0.9802 


0.9900 


0.9745 


O.9892 


0.9959 


0.9824 


0.9939 


0.9983 


0.9874 


0.9964 


0.9993 


0.9908 


0.9977 


0.9997 


0.9932 


0.9984 


0.9999 
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3 

3 

6 

6 

6 

9 

9 

9 

3 

3 

3 

6 

6 

6 

9 

9 

9 



TABLE 9 


PERIODIC 


-FLOW HEAT 


EXCHANGER 
C /C„<v, 


CO 

8 

Eh 

1 






r 




min 






1 




1 


2 


1 


2 


1 


2 


(hA)» 


^minA'max 


= 1.0 


Cmin/Cmax 


= 0.95 


Omin/Omax 


= 0.90 


1.00 


0.6672 


0.7280 


0.6760 


0.7406 


0.6846 


0.7530 


0.50 


0.6676 


0.7282 


0.6752 


0.7404 


0.6826 


0.7524 


0.25 


0.668U 


0.7284 


0.6751 


0.7404 


O.6815 


0.7521 


1.00 


0.7592 


0.8525 


O.7698 


0.8484 


0.7795 


O.8656 


0.50 


0.7597 


0.8525 


0.7687 


0.8481 


0.7769 


0.8628 


0.25 


0.7608 


0.8524 


0.7685 


0.8478 


0.7751 


0.8622 


1.00 


0.8017 


0.8765 


0.8129 


0.8940 


O.825O 


0.9099 


0.50 


0.8021 


O.8765 


0.8115 


O.8956 


0.8198 


0.9091 


0.25 


0.8052 


O.8765 


0.8m 


0.8933 


0.8176 


0.9084 




^mln/^max 


= 0.80 




= 0.70 


^mln/^max 


= O.50 


1.00 


0.7007 


0.7769 


0.7155 


0.7995 


0.7402 


0.8599 


0.50 


O.696O 


0.7752 


0.7075 


0.7964 


0.7245 


0.8551 


0.25 


0.6919 


0.7742 


0.7001 


0.79^3 


0.7099 


0.8270 


1.00 


0.7967 


0.8910 


0.8104 


0.9139 


0.8284 


0.9469 


0.50 


0.7902 


0.8888 


0.8000 


O.9105 


0.8102 


0.9398 


0.25 


0.7848 


0.8872 


0.7908 


0.9074 


0.7953 


0.9334 


1.00 


0.8594 


0.9365 


0.8514 


0.9560 


0.8659 


0.9789 


0.50 


O.8525 


0.9544 


0.8402 


0.9529 


0.8465 


0.9735 


0.25 


0.8265 

* 


0.9328 


O.8506 


0.9505 


0.8527 


0.9684 
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TABLE 9 CONTINUED 

C r / C min 



I 

1 



NTU 0 


(HA)' 


^min/^r 


5 


1.00 


0.7577 


5 


0.50 


0.7358 


5 


0.25 


0.7131 


6 


1.00 


O.836O 


6 


0.50 


O.813O 


6 


0.25 


0.7957 


9 


1.00 


0.8671 


9 


0.50 


0.8472 


9 


0.25 


0.8327 



2 


1 


2 


= O.30 


C min/ C max 


= 0.10 


0.8729 


0.7680 


0.8981 


0.8608 


0.7373 


0.8787 


0.81*84 


0.7136 


0.8579 


0.9658 


0.8379 


0.9745 


0.9546 


O.8133 


0.9597 


0.9439 


0.7957 


0.9457 


0.9881 


0.8675 


0.9908 


0.9808 


0.8472 


0.9821 


0.9734- 


0.8326 


0.9738 
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30 




31 




32 




33 




3 ^ 




55 




36 





37 




58 



10 



8 

6 

4 



AE/E„ 

(%) 



1 

0.8 

0.6 

0.4 



0.2 



0.1 



8 10 



Cr/Cndr 



FIG. 13 AE/Eoo vs Cr/Cmin In the range: 
1.0 ^ Cm-tn/Cmax ^ 0.90 
3.0 < NTO 0 ^ 9.0 
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9. Appendix 1 - Derivation of Iterative Equations 

Equations [2] , [3] , [ k] and [5J concerning the element on the 
side of C mflv , Fig. 2 (a), from the body of the thesis are repeated here 
for convenience. 



0 - " /xo) . • • 

Q = C, (Tft -Tri) (*M ‘ ' 

Q ^ ChA)„ AT^ ('M • • 

-- (WC* +T yo ) -C’/o)(r^ fT w ) 



, W 
. />7 

■ W 

• f5j 



If A^avg from equation [5] is substituted In equation [kj and the 
value for Q thus obtained is substituted in equation [2J the following 
will be obtained: 

(kX^Ur,; +i;. -7 rt )O r /w<) = - T, t )V/n,) 



This may be solved for *xo, giving 








%N X C wax, -"(KA)h 
X AAt £|>iQ*t (,WA')x 





In a similar manner T ro may be solved for by using equation ^3j/ » giving 



r f0 



_ I _ \(T- + T \ + I Xf/< Cr ~ I 77/ ■ ■ • ' 

~ Cr t(M). J 



At this point four relations can be derived which will be useful 
in subsequent manipulations to get expressions In terms of the dimen- 
sionless parameters. From the definition of NTU 0 it is obvious that: 
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I 



I 



■ • /'5 J 



^ nn>y\ 

ChA)n. 



NTTU 



n 



h/TUa [/+(hA)‘] 



An expression for reciprocal NTO X may be derived by multiplying and 
dividing it by NTU n an expressed in equation [15J , thus. 



1 1 C M)n, 

” a/tUx ~ Cm)* a/tm 0 c mi ^ 





CM)' 


["■] 




urujiH^y] c ^ : Vc^ A < 


The other two expressions, 


the derivations of which sure obvious. 


are 


as follows: 






<2. 




f'7j 




rJTUo [li M'J 


and. 


- 






Chfi) 


• [/?] 


II 

II 


tJTU 0 [/f M’J 


To obtain T xo in terms 


of only the inlet temperatures substitute 



equation [l 4 ] in equation £13] and solve for ^xo • Starting with equa- 
tions [13] and [lUj slightly rearranged, the algebra involved follows. 



\M±_ 

I NTH* 



+ / 









fo 



'll 

[Ami, 'J 



' X c 





I 



A A/* _ . * 



(T« + TJo) + 



Substituting [lba] in [l 3 a] > 



*v.t , I lr r ; ■■ ■[">*■] 

CWAV J 



N/Ttfx 



+ 



•Mt Cl + | ^X<J 

(hA)* J 




+ ”^£ 4 - 

(M)* 
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Each of the terms in the brackets may be put over the same common denom- 
inator which may then be multiplied out, giving 



haM'/.o .1 fliM. + iii _ -i+/]r„- 

*< 711 , (M), «ll, (M)« J r <» * WTO, CM), J 



Jjjj_ _ 


XA/*Cr 


ntu* 


(MX 


taV 


] rr ‘ 



Multiplying through by (hA)'/2N x C r and adding and subtracting 2T X ^ on 
the right gives 

A/r ,1 __ + N/y ,*] 1 ( t~ . t \ 

* 777; +l J r » - [«nu Cr V* T a ■ z( -7z-r r< -) 

Now solving for T x0 gives a fora convenient for iterative calculations, 
thus, 






t + 



Cimag hJf A/y 1 ^ 



Cr ^ A/7K 

The term in brackets is constant for a given problem so that 



"j Cr*/ ** I ri ^ > ['$] 







where K A may be put directly in terms of the dimensionless parameters 
giving 



< ‘-I / + 



Cm AX 



/ 

cxT 

Cvr,,‘rt 



tu. 

/Or 



fi/rtiJ/ + (hA)'] 





To obtain T ro in tferms of only the inlet temperatures substitute 

equation [13] in equation [lU] and solve for T rQ . Starting with 

equations [13] and flhj slightly rearranged, the algebra involved 
\ 

follows . 



X. 





z 



3 JJ_x 

h/TUn. 



i- I 



I XL 




kk 



ZfJ* Cc 
CfcA)* 




7~ro 



Ui + Ua 



+ 



Cr 

U/O* 




Tri 




Substituting [l3t>] in jl^bj , 



W* c «• 

(M)> 



+■ / - 






f I 



A/k 



Tr 0 ~ I (hA), 



"/ +■ 



I 



h7l4* 



T r ; 



+ 2 - 



i^r +i , 

WTUy J 



Ty/ 



Multiplying out the common denominator of all terms in brackets gives 

a&NfCr 4 A UCr Ml + i-/7r« - + TT - -I ill Tri 

( NTH* OiA)x A' 72 *' J ChA),VTKx ^ h/, ) x J 



+2-z) 7J/ 

NTH* J 



Multiplying through by NTU x /2N r and ending and subtracting 2T ri on the 
right gives 

‘ ^ 



Cr ^ Cr 
(h/0* N* Cnnccx 



•] 



+ ' T r 



A/x Cr- + 



ro 



( h A)> A/< 

t Z { lx i - / ri / 



Now solving for T rQ 
thus, 

I ro " I ri 



gives a form convenient for iterative calculations, 



+ l[ 



/ + 



Cr Wy 







M 



The term in brackets is constant for a given problem so that 



I <0 ~ I rt ^ 



t ( 2 (T% t ' — / ri ^ 



M 



^5 



where Kq may be put directly in terms of the dimensionless parameters 
giving 



4 Mi (m) ( 2 , m| 'y, 

m 0 [ i+M'J 



Cmin Gr ^ 



<* = *1 1 f + 



-i 



N 



For the element on the side of Cj^ , Fig 2(b), a set of equations 
similar to equations [2] through [5] can be written. They are: 



■mirt C"^l> "" ^ 1<’ /) C ) ' 


■ ■ ■ uo 


C< (,7V/ “ ^ro) (, - 


. . , • M 


ati^C 1 /^*- • • ' ‘ 


, . C 2 *] 


- CA)(Tr; +Tr«) “ {'/^(.Tni + 1 * 0 ) 


. . . Z' 2 ^ 



Following the same procedure as before results in two equations 
similar to equations [13J and [lUj as follows 

( bA) h 



T » f= l/r -r \ _ 

T "° - <»* J ( Trt + + [mu. V<m^ J t * ; 



■ Us] 



and 



T, 



( a)» 



+CM)v 



3lV»v d r - 
2/V* A- (M)n 



Tv 



V*L 



+ 1 ho) + 

These may be solved as before for the outlet temperatures in terms of 
only the inlet temperatures giving 



. l 24 >] 



r„„ = ui - u 1 + 



l 



T 7 ~ ^ 



.? A Ijr 

tfTK* 



-/ 



(Jni-T r i) . ... M 



lv6 



and 




/ ri 
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A/ f 



(.M)* _ 




The derivation of equations [27] and [ 28 J is also readily appar- 
ent by comparison of equation [25] with [l3j and [26J with [l 4 ] , The 
only difference is the substitution of the subscript "min" for "max" 
and "n" for "x" and then equation [27] follows from [19J and equation 
[ 28] follows from [20J . 

Equations [27 J and [ 28 J may be written as 



and 



T no - I ni + ^5 ( ( re " Ini J) , [ 7a J 



T r# = Tyi - K+ (T/i - T„i) ^ 



where 



K 3 - 2 



I 2AV -} 

' + Nr + W-mJftM)'] 



-1 
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/ ili- 



l Cr Mr. a. ~ ' '* '~ wsvv - 

K* = z [ 1 f 1 a/tu 0 [i+(m)] 



. , . . f»0 



. .... [f'l 



Equations [l 3 [j , [l 4 ] , [25J and [26 J are repeated here for con- 



venience. On the side of C, 



max ’ 



(M)* * 

,xo “ |2 V,c m ^ + (m) x 



U = 
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and on the side of n : 
"7" no ~ 

T"ro ~ 



Y 

t 


c5 A/ r Cmm 


ZtifC m ‘ lr , + (hA)^J 7~rt, J ^ 


t (M)*_ 



(.MX 1 _ v 1 JfitnC, - Mn 

+ (U)n 



khdia ■ 1 _ v 

av K e r + (hA)*J ''»<■ + y «^ + 



Tnc , . .[*5] 



7W , . . M 



It is obvious that the first term on the right of each of the fore- 
going equations is a positive quantity. It can next be observed that 
the coefficient of the second term of each equation must also be positive 
for if it were not the outlet temperature in question would grow hotter 
as its inlet temperature became colder, or vice versa, which would be 
contrary to the second law of thermodynamics. 

This provides a convenient set of conditions to be met before the 
iterative scheme can be expected to converge. By using the relations of 
equations [l5 ] through [l8] these conditions may be stated as follows: 



HTU.[it ‘AmyI 



Qrvuv\ 
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In the ranges of parameters considered and hy using an equal num- 
ber of subdivisions of each stream it can be seen that condition ^9_J is 
the most severe and therefore it alone can be used as the condition for 



convergence 



10. Appendix 2 - Computer Program Details 

The computer used for this work was the National Cash Register 102A 
general purpose digital computer and what follows assumes the reader has 
some familiarity with Its operation. 

The basic arrangement of the computer program was designed to accom- 
modate any number of subdivisions of each of the three streams up to and 
Including 64. It was found that, on the average, five iterations were 
required to meet the reversal condition and limit on the heat balance 
error with a Judicious estimate of the initial matrix temperature dis- 
tribution. With 16 subdivisions of each stream the time required for 
each iteration was four minutes. For different degrees of subdivision 
the time required for each iteration would vary, as would be expected, 
as the square of the factor by which the subdivision was changed, e.g. 

32 subdivisions of each stream would require 16 minutes for each iter- 
ation. It was determined that 16 subdivisions of each stream could be 
expected to give accuracy to about five units in the fourth place of the 
effectiveness. To finish in a reasonable time and obtain the desired 
four place accuracy it was decided that all the results would be cal- 
culated using 16 subdivisions of each stream to be followed by suffic- 
ient check calculations with 32 subdivisions to permit extrapolation of 
the results to infinite subdivision. 

The use of the various cells of the memory was as follows: 



Cells 


Use 


0000-0477 


Program 


0500-0577 


Initial T r *s 


0600-0777 


-Program 


1000-1077 


Current parameters and K's 




temporary storage. 
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1100-1177 



Outlet temperatures on side of Cma* 



1500-1377 



Outlet temperatures on side of 



1500-1577 



Initial-Intermediate-Final T r 's 



The composite computer program flow diagram is shown in Fig. 14 and 
the contents of cells 0000-0477 and 0600-0777 are recorded in Table 10. 
Box 10 and 11 of Fig. 14 accomplish the iteration process. This is 



of Cnjin in box 11. Since the internal temperatures are of no concern 
in this problem it was only neces saury to provide permanent storage for 
the fluid outlet temperatures of each column and the initial and final 
rotor temperatures of each row. 

Box 12 determines if the reversal condition is met by taking the 
difference between the initial and final T r of each row. A difference 
greater than the specified limit will cause the iteration cycle to 
start over. Also the largest difference is selected so that if the re- 
versal condition is met a means of manual inspection is provided, if 
desired, by causing this difference to be printed out upon putting switch 
2020 in the up position (box 17). 

If the reversal condition is met then the heat bala n ce error is 
computed in box l6 and If this is greater than the specified limit the 
iteration cycle is started over. The equation used to compute the 
heat balance error was arrived at as follows: 



done by the repetitive application of equations [ 6 ] to each element on 
the side of C mv in box 10 and equations [ 7 ] to each element on the side 



errof 
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where T no and T xo are the average fluid outlet temperatures . Since 
T xi vas taken as unity and T n ^ as zero then the heat balance error re- 



duces to 






Cr^ / _ \ 

St , I ' hrj J 



[ 30 ] 



(S'CfOf 



\ 



! (U o 



The average fluid outlet temperature on each side was computed by sum- 
ming the individual outlet temperatures of each column and dividing by 
the number of subdivisions. 

Another means of optional manual inspection was provided in box 19. 

If desired, the heat balance error would be converted from octal to decimal 
and printed out by putting switch 20^0 in the up position. 

By definition, the effectiveness may be computed by either of the 
following equations: 



where and T c0 are the average fluid outlet temperatures on the hot 
and cold sides, respectively. From the conditions assumed for calcula- 
tion purposes the following may be observed: 






[Silo] 




and therefore the expression for effectiveness from equation 




duces to 



E _ ( T no)avg 
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(T no ) a vg 11 have already been computed in the process of calculating 
the heat balance error and therefore in box 21 it is only necessary to 
convert this to decimal and print it out. 

The final T r 's punched on cards in box 22 were used as estimates 
for the initial T r ’s in subsequent calculations. This was usually done 
by using the final T r 's obtained at a given C^p/CmA* as the initial 
T r ' s for the next lower value of Cmin/Cmax* at corresponding values of 
the other parameters. . 

The remainder of the program consists of a provision for interupting 
the calculations at logical stopping places and a system of inserting 
new values of parameters sifter each solution. These are self-evident 
in Fig. 14. Parameters are entered as floating point numbers and the 
arithmetic involved in calculating the K’s is done in floating point. 

The K's, scaled down by a factor of two, sure then converted to single 
precision and the subsequent iteration is accomplished in fixed point, 
single precision. 

The basic computer program which is recorded in Table 10 is set 
up to use l 6 subdivisions of each stream. To modify it for 8 or 32 
subdivisions the following cells should be changed as indicated: 



Cell 


8 Subdivisions 


32 Subdivisions 


0307 


21050021000010 


21050021000040 


0^03 


5 


7 


0406 


35051021001510 


35054021001540 


0410 


36200015102001 


36200015402001 


0411 


35200021001110 


35200021001140 


0424 


36150720002001 


36153720002001 


0425 


36150720031507 


36153720031537 
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0430 


35200021001310 


35200021001340 


0433 


2ff0003 


2ff0005 


0436 


30131004332007 


30134004332007 


0437 


30111004332007 


30114004332007 


0443 


35151021000510 


35154021000540 


0447 


36051015102000 


36054015402000 


0456 


04300030000510 


04300030000540 


0462 


05300030001510 


05300030001540 


In its 


present form the program by-passes box 



to box 27. This causes all calculations to be made with (hA) ' = 1. To 



get results with (hA.)' = 1, l/2, and l/4 enter the following changes: 
0007 21100604020002 

0367 33100606760371 

0370 34300021000373 

0371 36100604001006 

0372 34300021000006 

In its present form the program produces results with C r /C mln = 

1, 1.25, 1.5, 2, 3, and 5. To eliminate Cr/C^n = 5 enter the following 



changes : 

0375 34067210020600 

0376 34067310030600 



Space is provided in the iteration subroutines for entering appro- 
priate commands to get any desired temperatures printed out before they 
sure lost. These may be printed in decimal simply by planing in cell 
2000 and unconditioned transfer to the decimal conversion and print sub- 
routine at cell 0700. 



34 




Fig. lb COMPUTER PROGRAM FLOW DIAGRAM 
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TABLE 10 COMPUTER PROGRAM COMMANDS 



0000 31014-0004011000 

0001 31040004011002 

0002 31040004011004 

0003 51040004011006 

0004 21100004020002 

0005 21100204020002 

0006 21100404020002 

0007 34300021000010 

0010 31040004012000 

0011 31100610072002 

0012 34300021000733 

0013 31200420052000 

0014 31100410052002 

0015 34300021000725 

0016 31200420052002 

0017 51040304012000 

0020 34500021000760 

0021 51200420051040 

0022 51200420052000 
0025 51100610072002 

0024 54500021000725 

0025 51200420051042 

0026 51100210052002 
0027 54500021000725 

0050 51200420051044 

0051 51040004012000 



0052 54500021000760 
0055 51200420051046 

0054 51104210452000 

0055 54500021000725 

0056 51200420051050 

0057 5H00010012002 

0040 54500021000760 

0041 51200420051052 

0042 51100210052000 
0045 54500021000725 

0044 51200420051054 

0045 51200420052002 

0046 51040004012000 

0047 54500021000760 

0050 51200420051056 

0051 51105210552000 

0052 51105610572002 
OO55 54500021000753 

0054 51200420051060 

0055 51200420052002 

0056 51040004012000 

0057 54500021000753 

0060 51200420052002 

0061 54500021000760 

0062 51200420051010 
OO65 51105010512000 
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TABLE 10 CONTINUED 



0064 31105410552002 

0065 34300021000753 

0066 31200420052002 

0067 31040004012000 

0070 34300021000733 

0071 31200420052002 

0072 34300021000760 

0073 31200420051012 

0074 31104010412000 

0075 31104610472002 

0076 34300021000733 

0077 31200420052002 

0100 31040004012000 

0101 34300021000733 

0102 31200420052002 

0103 34300021000760 

0104 31200420051014 

0105 31104410452000 

0106 31100210032002 
0107 34300021000733 

0110 31200420052002 

0111 31040004012000 

0112 34300021000733 

0113 31200420052002 
on4 34300021000760 
0115 31200420051016 



0116 30101110101011 

0117 30101310121013 

0120 30101510141015 

0121 30101710161017 

0122 33210010110010 
0123 33210010130010 

0124 33210010150010 

0125 33210010170010 

0126 17210030000640 

0127 34300021000130 

0130 35040421000131 

0131 35052021001520 

0132 35013104050131 

0133 34040601310131 

0134 35041221000157 

0135 35041321000141 

0136 35041421000146 

0137 34300021000140 

0140 35040121002000 

0141 36200015202001 

0142 30200104002001 

0143 252OOIIOJLL2OO2 

0144 25200110132003 

0145 36200020022000 

0146 35152020031520 

0147 34300021000154 
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TABLE 10 CONTINUED 



0150 00000000000000 

0151 00000000000000 

0152 00000000000000 

0153 00000000000000 

0154 35014104070141 

0155 35014604050146 

0156 34041001410141 

0157 35200021001120 

0160 34300021000170 

0161 00000000000000 

0162 00000000000000 
0165 00000000000000 

0164 00000000000000 

0165 00000000000000 

01 66 00000000000000 
0167 00000000000000 

0170 35015704000157 

0171 34041101570135 

0172 34300021000200 

0173 00000000000000 

* 

0174 00000000000000 

0175 00000000000000 

0176 00000000000000 

0177 00000000000000 

0200 35042721000221 

0201 35042421000204 



0202 35042521000211 

0203 35210021002000 

0204 36147720002001 

0205 30200104002001 

0206 25200110152002 

0207 25200110172003 
0210 35200020022000 

0211 36147720031477 

0212 34300021000216 

0213 00000000000000 

0214 00000000000000 

0215 00000000000000 

0216 36020404150204 

0217 36021104050211 

0220 34020404260204 

0221 35200021001320 

0222 34300021000231 

0223 00000000000000 

0224 00000000000000 

0225 00000000000000 

0226 00000000000000 

0227 00000000000000 

0230 00000000000000 

0231 35022104000221 

0232 34043002210201 

0233 34300021000314 
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TABLE 10 CONTINUED 



0254 00000000000000 

0235 00000000000000 

0236 00000000000000 

0237 00000000000000 

0240 35043421000243 

0241 35043521000255 

0242 35210021002000 

0243 30132004332007 

0244 35200020072000 

0245 35024304150243 

0246 34043602430243 

0247 30200004001040 
0250 31210020002000 
0251 31100010012002 

0252 34300021000725 

0253 31200420052000 

0254 35210021002003 

0255 30112004332007 

0256 35200320072003 

0257 35025504150255 

0260 34043702550255 

0261 36040120032003 

0262 31210020032002 

0263 35200204002002 

0264 34300021000265 

0265 34300021000760 



0266 30200520042005 

0267 36040120052005 

0270 30200504002005 

0271 34200504400300 

0272 34300021000326 

0273 00000000000000 

0274 00000000000000 

0275 00000000000000 

0276 00000000000000 

0277 00000000000000 

0300 35044221000300 

0301 35150021000500 

0302 35030104050301 

0303 34044303010301 

0304 17201030000306 

0305 34300021000134 

0306 21044404450001 

0307 21050021000020 

0310 22000000000000 

0311 34300021000134 

0312 00000000000000 

0313 00000000000000 

0314 35210021001050 

0315 35044621000316 

0316 36052015202000 

0317 34200010500324 
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TABLE 10 CONTINUED 



0320 35031604500316 

0321 34044703160316 

0322 34105004410300 

0323 34300021000240 

0324 35200021001050 

0325 34300021000320 

0326 17202030000351 

0327 17204030000335 
0550 34300021000340 
0331 21105021000001 

0332 34300021000327 

0333 35200521002000 
0354 34300021000700 

0335 34300021000340 

0336 00000000000000 

0337 00000000000000 

0340 34300021000347 

0341 00000000000000 

0342 00000000000000 

0343 00000000000000 

0344 00000000000000 

0345 00000000000000 

0346 00000000000000 

0347 35104021002000 

0350 34300021000700 

0351 21044404450001 



0352 35045721000355 
0555 05300030001520 

0354 13046004550000 

0355 30300004611077 

0356 13046004550004 

0357 30300004611077 

0360 35035504530355 

0361 34046203550355 
0562 34500021000363 

0363 34300021000364 

0364 17201030000366 

0365 34300021000367 

03 66 22000000000000 

0367 34300021000373 
0370 00000000000000 

0371 00000000000000 

0372 00000000000000 
0375 34065210040603 

0374 34065310050603 

0375 34067410020600 

0376 34067510030600 

0377 34300021000611 

0400 00000000000000 

0401 00400000000000 

0402 02000000000000 

0403 00000000000006 
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TABLE 10 CONTINUED 



0404 35050021001500 

0405 00000100000001 

0406 35052021001520 

0407 00000000010000 

0410 36200015202001 

0411 35200021001120 

0412 35200021001100 

0413 36200015002001 

0414 35150020031500 

0415 00000100000000 

0416 30050004002000 

0417 30052004002000 

0420 30152004002000 

\ 

0421 50150004002000 

0422 30110004002000 

0423 30112004002000 

0424 36151720002001 

0425 36151720031517 

0426 36147720002001 

0427 35200021001300 

0430 35200021001320 

0431 30130004002000 

0432 30132004002000 

0433 02000000000004 

0434 30130004332007 

0435 30110004332007 



6436 30132004332007 
0437 30112004332007 

0440 00000100000000 

0441 00000004000000 

0442 35150021000500 

0443 35152021000520 

0444 00341212341212 

0445 10000000000000 

0446 36050015002000 

0447 36052015202000 

0450 00000100010000 

0451 04300030000500 

0452 15555555555556 

0453 00000000000010 

0454 00000000000500 

0455 11111111111111 

0456 04300030000520 

0457 05300030001500 

0460 26666666666666 

0461 00000000001000 

0462 05300030001520 

0463 00000200020000 

0464 31066406651002 

0465 31210006551000 

0466 00000000000000 

0467 00000000000000 
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TABLE 10 CONTINUED 



0470 00000000000000 

0471 00000000000000 

0472 00000000000000 
0475 00000000000000 

0474 00000000000000 

0475 00000000000000 

0476 00000000000000 

0477 00000000000000 

0600 31066406651002 

0601 35060004650600 

0602 34300021000002 

0603 31040004012000 

0604 31100410052002 

0605 34300021000733 

0606 31200420051004 

0607 31040004011006 

0610 34300021000006 

0611 35046421000600 

0612 31100010012000 

1 

0613 35210021002002 

0614 36210OO6632OO3 

0615 34300021000735 

0616 34200521000621 

0617 35046521000621 
0620 22000000000000 
0621 31210006551000 



0622 35062104070621 

0623 34300021000001 

0624 00000000000000 

0625 00000000000000 

0626 00000000000000 

0627 00000000000000 

0630 00000000000000 

0631 00000000000000 

0632 00000000000000 

0633 00000000000000 

0634 00000000000000 

0635 00000000000000 

0636 00000000000000 

0637 00000000000000 

0640 35045121000646 

0641 05300030002100 

0642 06045204550000 

0643 30300004541077 

0644 06045204550004 

0645 30300004541077 

0646 04300030000520 

0647 35064604530646 

0650 34045606460641 

0651 34300021000130 

0652 00000000000004 

0653 00500000000000 
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TABLE 10 CONTINUED 



0654 00000000000001 

0655 007^31^631^63 

0656 007146314631^6 

0657 00631463146315 

0660 00546314631463 

0661 00400000000000 

0662 00231463146315 

0663 00063146314632 

0664 00000000000001 

0665 00500000000000 

0666 00000000000001 

0667 00600000000000 

0670 00000000000002 

0671 00400000000000 

0672 00000000000002 

0673 00600000000000 

0674 00000000000003 

0675 00500000000000 

0676 02000000000001 

0677 00400000000000 

0700 27300007162001 

0701 35071420010713 

0702 35200021000722 

0703 05300030000714 

0704 32200207002000 

0705 25200220072002 



0706 32200220012000 

0707 32210020012002 

0710 27200120052001 

0711 34200121000705 

0712 21200007200001 

0713 34300021000351 

0714 34300021000000 

0715 02000000000004 

0716 02000000000030 

0717 00500000000000 

0720 01000000000000 

0721 00740000000000 

0722 00463511721510 

0723 02000000000001 

0724 00000000000001 

0725 30300007162004 

0726 35071420040732 

0727 35200020022004 

0730 25200120032005 

0731 31200420052004 

0732 34300021000253 

0733 30300007162004 

0734 35071420040745 

0735 33200020020746 

0736 36200020022004 

0737 30200120042005 



63 



TABLE 10 CONTINUED 



07^0 35200321002006 

0741 35200221002004 

0742 35200520062005 

0743 37200521000755 

0744 31200420052004 

0745 34300021000113 

0746 36200220002004 

0747 30200320042005 

0750 35200121002006 

0751 35200021002004 

0752 34300021000742 
0755 27200507232005 

0754 30200507232005 

0755 27200507242005 

0756 35200407242004 

0757 34300021000744 

0760 30300007162004 

0761 35071420040766 

0762 36200020022004 

0763 23200120032005 

0764 37200521000767 

0765 31200420052004 

0766 34300021000266 

0767 27200507232005 

0770 30200507232005 

0771 27200507242005 



0772 35200407242004 

0773 34300021000765 

0774 00000000000000 

0775 00000000000000 

0776 00000000000000 

0777 00000000000000 
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